A mode-coupling theory for the pasty rheology of soft glassy materials 
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In this paper, we introduce a simple mode-coupling model for concentrated suspensions under 
flow. This model exhibits a jamming transition, and stress vs shear rate relations which are very 
similar to experimental results. Namely newtonian regime or yield stress are followed by a slow 
variation of the stress for higher shear rates, and by an apparent newtonian regime for very large 
shear rates. Another striking result is that under oscillating strain, even in the jammed state, the 
system exhibits a relaxation time which depends on the strain amplitude. 
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PACS numbers : 62.20.Fe, 83.20.-d, 83.50.Gd 



We present a model which aims at describing the mechanical behavior of very concentrated suspensions of soft 
i—i particles. Such systems, like concentrated emulsions, colloidal suspensions near random close packing, concentrated 
suspensions of small particles of gel, and concentrated colloids electrostatically charged, are known to exhibit peculiar 
paste-like behavior for high enough concentration of the dispersed phase : they have a yield stress and thus very 
non-linear mechanical properties. They all share the property of being made of a collection of soft objects. 
Beyond a critical value of their concentration, the particules become packed together, and they experience a jamming 
transition. These systems have recently been named Soft Glassy Materials. 

At the macroscopic scale, the jamming transition corresponds to the fact that, below a critical concentration, the 
system is fluid, and the stress vanishes at rest ; but, above the critical concentration, stress can be stored in the 
system — after flow cessation for instance — . Obviously, because these systems exhibit a yield stress, their mechani- 
cal behavior under the action of flow is very non-linear. More precisely, some recent model systems and simulations 
exhibit a power-law relation between stress and shear rate, with quite small exponents — around .1 — for small shear 
I ' rates, and an exponent around 1 for large shear rates 
, At the mesoscopic scale, the jamming transition means that the motions of the particles under flow become more 
and more collective as it is approached, leading to very complex trajectories and rearrangements of the particles. 

t>; 

' Actually, this jamming situation has two generic ingredients : 

o : 
oo : 

• when a flow is applied to the system, it modifies deeply its dynamic. 

> ; 

The effect of flow on a system near jamming transition may be described by a phase diagram shear rate vs con- 
Q centration. At zero shear rate, random close packing concentration, <fi c , corresponds to a critical point. Below (j) c , 
any infinitely small applied stress induces rearrangements and flow of the particules, whereas above C , a finite non- 
vanishing stress must be applied in order to make the system flow. This means that there is a half-line of degeneracy 
of the stress — for vanishing shear rate, above the critical concentration — , with a critical point at its extremity. But 
O ■ as far as the shear rate is different from zero, the stress always reaches a unique value in a stationary flow, whatever 
, the concentration is. In other words, the shear kills the jammed state, because the shear rate renews constantly the 
structure of the system . This description remains valid as long as the stress remains homogeneous in the sample 
(no fractures for instance) . This phase diagram is very similar to the magnetic field vs temperature phase diagram of 
the Ising model. The shear rate 7 plays the role of the magnetic field, the stress a that of the magnetization, and the 
concentration cp stands for the temperature. However, Soft Glassy Materials are not at thermodynamic equilibrium, 
and the relations between 7, a and (f> do not derive from a free energy, but from a dynamical equation. Here, we 
introduce such an equation, based on a mean-field analysis of a mode-coupling theory. 



at rest, above a critical concentration, the system is jammed ; 
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An attempt of this type of approach was recently proposed by P. Sollich, M.E. Cates and us M. The jamming 
transition was described using a glass model introduced by J. -Ph. Bouchaud ||. One of the results of the model was 
the scaling law relating the stress with the shear rate, with an exponent depending on the distance to the jamming 
transition. The control parameter was an effective temperature which attempted to describe the mechanical noise. 
It was introduced phenomenologically, and was not related to the macroscopic shear rate, whereas shear must induce 
changes of configurations in the sample, and thus affect its effective temperature. 

In this letter, we present a model which describes more directly the interactions between the particles trajectories. 
This model thus belongs to the mode-coupling class of models B . 

We first divide the sample into blocks carrying a stress a depending both on time t and on their position r. The 
evolution of the blocks follows simple rules : 

1. When submitted to a shear rate 7, the stresses of all the blocks increase during time interval dt by the quantity 
Gajdt. 



2. If the absolute value of the stress of a given block is larger than a critical value er c , it is set to zero after a time 
t. This processus is very similar to rearrangements described by Princen for foams @]. However, the fact that 
the stress is set to zero after the rearrangement, or randomly distributed around zero is clearly not critical. 



3. If the stress on r' is set to zero, stress field is modified in the whole sample. We will describe this interaction 
with a mean field approximation, by adding a random value A to the stresses of all the other blocks . 

Let us discuss the last point. Actually, if the block in r' relaxes from a to zero, the stress at point r varies of the 
quantity Q(r' — r)er(r) where Q is the stress/stress elastic propagator. Q is well known in an homogeneous system 
and it appears that its mean value for the 4-7T steradians is equal to zero. This implies that the random value A must 
be chosen with a mean value equal to zero. In an homogeneous elastic medium, the stress propagator is short range 
- it decreases as 1/r 3 — , and very anisotropic, leading to step by step propagation of the stress along preferential 
directions, and finally fracturation phenomena, as described in earthquakes models 

But, systems we are interested in are not at all homogeneous, so the stress propagates along easy paths defined by 
successive contacts, and stress thus propagates over a long range. Moreover, these stress paths are very fragile and 
renew completely after very small variations of the strain, so mean field approach may be a good way to describe 
these systems. 

Finally, we can write the evolution for the probability P(o~, t) of finding a stress a in a block, at time t : 

d t P(a, t) = -GojdrPia, t) + Dd 2 a2 P(a, t) - g(H ~ ^ P(a, t) + - f P(a, t)daS(a) (1) 

T T J\<j\>a c 

where H(x) = 1 for x > 0, and otherwise, and 5 is the Dirac function. 

The first term corresponds to the first rule, the second one to the noise of the third rule, and the last two terms 
express the failure for stresses above a c with a rate 1/r and the reset of the stresses to zero after the failure. The 
third condition writes : 



D = a\\ P(<T,t)da (2) 



r 



where D is the amplitude of the noise, and is then proportional to the density of blocks that rearranged during time 
t, by a proportional factor a, that depends a priori on the microscopic properties of the material. It could represent 
something like a mechanical fragility. 

Let us first describe the solution of this model in the absence of shear. We get a naive description of a mode-coupling 
theory of glass transition |Q. The model exhibits a dynamic transition, between a frozen system, and a liquid state 
in which the stress is self-sustained, through a diffusive term, by over-stressed regions. This description of a jamming 
transition is very similar to the one due to Ivanov et al. flfj|| . 

For 7 = 0, one easily finds that (||) and (g) leads to : 

„ cuD , . 

D = = 3 

1 + 2VD + 2D 
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For a < | = a c , the only solution is D = 0. Solutions correspond to distributions of stable states of the local stress, 
i.e. distributions in which P is equal to zero for |er| > a c . This is quite obvious : all the terms in (|l]) vanish in that 
case. 

On the opposite, for a > a c in addition to the previous solution, it exists a stationary solution with a non-zero 
value of D. The system can continuously rearrange as a function of time. It is easy to see that D scales simply as 
(a — a c ) 2 above a c . 

Let us now describe the solution in the presence of a stationary shear rate, 7 = c te . In that case, equations (□) and 
have always a single stationary solution (Fig. [l]). The derivation is straightforward. We first solve equation (jl]) for 
a given value of D, and then write the self-consistent equation for D i.e. eq. (|^). Then, we make expansions of this 
self-consistent equation to get the scaling laws given further. From the distribution probability P, we compute the 
macroscopic stress, cr macrol via the mean field assumption : the macroscopic stress is approximated by the average 
value of the local stresses. 

Whatever the value of a, at high shear rates (7 > 1/r), the flow is Newtonian, and the viscosity scales as the 
life-time r of the excited state multiplied by the elastic modulus G . Collective effects are not important in this 
situation. 

For small values of the shear rate, the system exhibits a Newtonian regime in the liquid phase, and a yield stress in 
the jammed phase. The yield stress appears because the solutions of (|l|) and (Q) are degenerate for a vanishing value 
of the shear rate, leading to a discontinuity of the shear stress as a function of the shear rate. 

We find the following exponents, in the limit of very small shear rates : 

for a > a c , a ~ G A fr{a — a c )~ 2 (4) 
7— >o 

for a < a c ,a ~ G(a c - a) 1/2 (5) 
7^0 

For the critical value a = a c , we find a scaling between the stress and the shear rate: 

for a = a c , a ~ 7^ 5 - (6) 

7^0 

Moreover as far as a is in the vicinity of a c , the stress always scales like 7 1 / 5 in the crossover regime between the 
high shear rate and low shear rate limits. (Fig. |[ ||). 

These many regimes mimic the recent experimental results of M. Cloitre on microgel beads around the close packing 
of the beads, and that of Mason et al. on concentrated emulsions [jll). It describes also the results of D. Durian's 
simulations of foams above the close packing, and some of his recent experiments. The value 1/5 of the exponent 
may seem odd ; it's not so far from the experimental results of Cloitre (0.1), and from simulations of foams using the 
Durian model by A.Liu (0.16) |L2j, as well as the experimental results on foams by Durian (.15) 

The main other rheological properties of concentrated systems may be studied thanks to the response to a periodic 
flow, 7(i) = 70 cosui, expressed in terms of the elastic and loss modulus. In order to calculate the elastic and the loss 
modulus for a non-vanishing amplitude of the strain, we make an expansion of the evolution equation to the second 
order in 7. The calculation is analytical, and we only give the final results. 

First, we obtain that the response is Maxwellian at low frequency, whatever the values of a and 70. In other words, 
G' and G" scale respectively as 10 and uj 2 when uj — > 0. 

At high frequency, the behavior is more complex : whatever the phase, G' is constant, whereas G" exhibits a 
maximum, uj c . It decreases as w -1 / 2 in the liquid phase, and as uj~ 2 in the glassy phase. The more striking result 
lies in the dependance of the characteristic frequency oj c with a, which depends on the amplitude of the deformation. 
Indeed, for vanishing amplitudes, uj c — in the glassy phase, and uj c ~ (a — ot c ) 2 in the liquid phase. This is quite 
obvious, as the high frequency modulus is nearly constant, uj c reflects the behavior if the viscosity. But, for non-zero 

amplitudes, u) c ~ I7I in the glassy phase, and oj c ~ (a — a c ) 2 + 7 in the liquid phase (Fig. 4). 

7— >o 7— »o 

So, in the glassy phase, a periodic deformation induces a continuous flow of reorganizations playing the role of an 
effective temperature. This induces a relaxation process with a characteristic time scaling as |7| _1 - It thus exists in 
this phase infinite non-linearities, with an apparent slow relaxation. 

The important point lies in the result that a finite strain amplitude modifies completely the divergence of the typical 
time of mechanical relaxation. This has to be compared with experimental results showing that pasty systems exhibit 
very non-linear viscoelastic behavior . Our model leads to an infinitely small linear regime in the jammed state. 

So, our model describes in a very naive way a jamming transition, and mimics many of the experimental results on 
soft glassy materials. We hope that this primitive model opens new gates to further more subtle descriptions of the 
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jamming transition. Questions arise whether there is some universality or not in these systems. The experimental 
situation is not clear, as most of the mechanical measurements do not deal with the jamming transition itself and its 
scaling. Moreover, the spatial homogeneity of the flow depends dramatically on the systems : some of them fracture 
macroscopically in flow, while other do not, for completely unknown reasons , and mechanical properties are very 
inhomogeneous. More than ever, precise spatial averagings of the mechanical properties are really a challenge in these 
jamming systems under flow. 
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FIG. 1. Distribution probability P(a) of the stress, for different values of a. 



FIG. 2. Stress a vs shear rate 7 for different values of a. 



FIG. 3. Dynamical phase diagram : yield stress a c vs a for different values of the shear rate 7. 



FIG. 4. Behavior of the characteristic frequency, lo c vs a at constant peak amplitudes of the shear. 
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